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A NUMERICAL METHOD FOR EVALUATION AND UTILIZATION 
OF SUPERSONIC NACELLE-WING INTERFERENCE 


By Robert J. Mack 
Langley Research Center 

SUMMARY 

A numerical method for calculating (1) the inviscid flow nacelle-wing interference 
loads on a warped wing surface, (2) a load- compensating reflex camber surface for opti- 
mization purposes, and (3) the interference -lift effective -area contribution to sonic -boom 
pressure disturbances has been presented and discussed. Within the limitations of a 
modified supersonic linearized theory, the numerical method can provide reasonable good 
estimations of the loads imposed on a wing surface by nacelles mounted nearby. 

INTRODUCTION 

Application of supersonic flow theory to the design and evaluation of aerodynamically 
efficient supersonic-cruise aircraft has been greatly facilitated by the recent employment 
of the high-speed digital computer. Simple numerical methods which closely approximate 
the behavior of complex integral equations can be programed for the computer to give 
quick solutions to aerodynamic theory problems and can extend greatly the range of prob- 
lems that may be handled. For example, the program described in reference 1 computes 
an optimized camber surface at a given Mach number for a wing of arbitrary planform, 
whereas those of references 2 and 3 determine the aerodynamic characteristics of a given 
wing with an arbitrary planform and camber surface. The wave drag of an entire con- 
figuration may be evaluated with the program outlined in reference 4. Similar unpublished 
programs based on the Sommer-Short T-prime method permit the calculation of the skin- 
friction drag of the configuration. With these programs, the overall aerodynamic char- 
acteristics of a proposed design can be found with a reasonably good degree of accuracy 
for most slender -body thin -wing aircraft. 

In estimating the overall aerodynamic characteristics of a configuration employing 
nacelles or stores, a means of evaluating the interference forces is required. The com- 
putation technique described in this report provides a method for determining the distri- 
bution and intensity of the interference forces on a warped wing surface. 

However, airplane design studies which must include nacelle -induced drag minimiza- 
tion require more than a knowledge of the interference forces. Achievement of the highest 
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possible level of aerodynamic efficiency demands careful consideration of the manner 
in which the nacelles and the wing are integrated. Current design practice, as discussed 
in reference 5, is based on the concepts of favorable interference and wing reflexing 
(warping the wing surface away from the nacelles). Nacelles positioned under the wing 
trailing edge give, in general, relatively low wave drag. Wing reflexing utilizes the idea 
that a localized warping of the wing surface will produce an incremental loading opposite 
of that induced by the nacelles and will give a combined lift and moment loading, due to 
wing and nacelle, which can be the same as the optimized loading originally designed for 
the wing. A numerical method for determining reflexed wing coordinates which makes use 
of the previously discussed nacelle interference loadings and incorporates the wing cam- 
ber design method of reference 1 is considered in this report. 

Also discussed is the effect of the nacelle-wing interference on the flow field under 
the aircraft. Changes in the airplane lift distribution brought about by the nacelle pres- 
sure field acting on the wing surface affect not only the drag, but also the distribution of 
the pressure in the shock field surrounding the airplane. Thus, the nacelle flow-field 
computation is also useful in a sonic -boom analysis using the methods discussed in 
reference 6. 

In summary then, these three items are discussed in this report: (1) calculation of 
the nacelle-wing interference forces, (2) calculation of a reflexed wing camber surface, 
and (3) calculation of the incremental lift distribution which is used in a sonic -boom 
analysis. Necessary assumptions are mentioned as they are needed in the discussion 
and presentation. 

Since the object of the present study is to provide a simplified method of accounting 
for nacelle effects to supplement existing configuration analysis programs, the use of more 
refined or more rigorous solutions than those employed in the previously mentioned work 
would seem to be unnecessary. However, the calculation of interference forces and wing 
reflexing requires that very good approximations to actual pressure disturbance locations 
be made. The method used in this report utilizes the modified linearized theory developed 
by G. B. Whitham (ref. 7). It not only overcomes the difficulty of obtaining solutions for 
ducted nacelles, but also calculates the positions and strengths of shock waves originating 
at the nose and surface discontinuities of the nacelle or store. Several examples are 
presented to show comparisons of nacelle pressure fields calculated by the numerical 
method with those measured in wind-tunnel experiments, and to illustrate practical 
applications to aircraft design analysis. 
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SYMBOLS 


A nacelle cross-section area 

A(L,N) interference zone grid element weight factor 

b wing span 

B(L,N) wing area grid element weight factor 

c wing chord 

c mean geometric chord 

Cj) pressure drag coefficient 

Cl lift coefficient 

C m moment coefficient about 0.25c based on mean geometric chord 

C m>0 moment coefficient at Cl = 0 

Cp pressure coefficient 

ACp lifting pressure coefficient 

F(t) area distribution function (called F(y) in ref. 7) 

h area influence function (called h(x) in ref. 7) 

K reflex factor, Load canceled 

Load induced 

l maximum length on wing 

l n nacelle length 

L,N designation of influencing grid elements (see fig. 3) 

L*,N* designation of field-point grid elements (see fig. 3) 
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M 


Mach number 


x,y,z distances along Cartesian coordinate axes X, Y, and Z 

x,r distances along cylindrical coordinate axes X,R 

t dummy variable along X-axis, used in F(t) calculation 

R body radius 

S wing area 

Uoq free- stream velocity 

P = \/M 2 - 1 

y 

A 

<P 

Subscripts: 
av average 

c camber surface 

e effective 

I due to interference 

max maximum 

min minimum 

si shock-wave — camber -surface intersection 


ratio of specific heats, 1.4 for air 
dummy variables along x,y directions 
leading-edge sweep angle 
perturbation potential 
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le 

leading edge 

te 

trailing edge 

n 

value on nacelle 


Primes denote derivatives with respect to length t. The symbol A denotes 
incremental quantities. 

THEORY AND DISCUSSION 
Interference Force Determination 

A simplified representation of the nacelle interference problem is shown in figure 1. 
The sketch depicts a warped-delta -wing ducted-nacelle configuration in supersonic flow. 

A fuselage is not shown so that attention can be focused on the nacelle-wing interference. 
The determination of interference forces on this configuration is made with due regard to 
the capabilities of other computer programs with which it is to be used concurrently in 
analysis work. Wave drag due to wing and nacelle thickness (including mutual interfer- 
ence) is accounted for in the wave drag program of reference 4. Therefore, the method 
of this report is applicable only in the calculation of the lifting forces on the wing resulting 
from nacelle thickness, and to any drag component that results from the inclination of the 
wing mean camber surface. The effect on the nacelle of the wing lift-induced pressure 
field has been neglected since that drag contribution is assumed to be small. Similarly, 
the effect of the reflection of the nacelle pressure field from the wing surface has been 
neglected. It might be noted that a wing reflex surface designed to compensate for the 
nacelle-induced wing loading will reduce or eliminate this latter drag contribution. 

Figure 2 shows the nacelle referencing coordinates, and figure 3, the grid element 
system used to define the wing and its camber surface. This coordinate system is the 
same as that used to describe the wing in references 1 to 3. 

The linearized flow equation for the perturbation potential about a body of revolu- 
tion in supersonic flow is 


( 1 ) 

9x z r 0r z 

Bodies with noncircular cross sections can also be analyzed if the geometry and probable 
attendant flow field will not be greatly distorted by using one or more bodies of revolution 
whose total area and area buildup are the same as those of the original (ref. 4). 
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The classical solution is derived by representing the body by a line of sources and 
having the pressure disturbances propagate outward along Mach lines. G. B. Whitham 
modified this solution by replacing the Mach line characteristic 

x - /3r = Constant (2) 


with a more exact characteristic 


where 


x - /3r = t - — + 1)M V F(t) 

1/2? 


(3) 


F ( t) = r° 


/3R(I) 


h 


dA\Z) 

pm)) 2 tt 


(4) 


and t is a constant equal to x - /3R(x) where the characteristic meets the body surface 
(eq. (12), ref. 7). This Stieltjes integral form for F(t) is used instead of the "smooth" 
body form 


F (t) = ^- 

2jt Jn 


A"(€)d4 


j 0 ]ft 


(5) 


because, in general, nacelles do not always have smooth and continuous surface slopes, 
and the effects of these discontinuities are more noticeable in the flow field near the body 
which impinges on the wing surface. In the case where the body shape is smooth, and 
surface slopes are analytically continuous, both forms of F(t) give the same solution. 

Ahead of the inlet, the flow is undisturbed by the nacelle; thus a long cylindrical 
stream tube can be extended forward from the intake without changing the flow character- 
istics. Along this stream tube, AA' = 0; therefore, F(t) = 0 for all t < -i3R(0). This 
initial point t = -/3R(0) is the intersection on the body axis of a Mach line from the 
nacelle lip and corresponds to the point where an abrupt change in source strength would 
cause a discontinuity in AA' on the stream-tube — nacelle surface. The step change in 
F(t) due to this flow deflection at the lip is 


AF(t) = 1 /— 2— (( 

|//3R(0) 2n 

■ which is equation (26) of reference 7. Beyond the lip source point, t > -/3R(0), values of 
F(t) are calculated from equation (4) in the form 
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( 7 ) 


where a summation operation replaces an integration over a range of | starting at 
£ = -/3R(0) and extending to the point on the body surface where a Mach line from t 
intersects (see fig. 4(a)). The field point distance r and the body radius R(t) deter- 
mine the slope AF/At of the lines which change the F(t) curve to a pressure profile 
through the relationship 


AF _ r + 1 M 4 (g) 

^ ^2/3 3 yfr - ^R(t) 

as seen in figure 4(b) and described in reference 8. Figures 4(a) and 4(b) show the com- 
ponents which make up the F(t) calculation, their relationship to nacelle thickness, 
nacelle surface slope and Mach number, and that a ACpj can be found by suitably 
modifying the F(t) curve. Notice that the approximation 


ft \/r - v® 

is very good when r is more than one or two body lengths, and the nacelle or pod is 
truly "slender." 

The calculation of shocks and shock locations from the F(t) curve is explained 
fully in reference 7, and the computational technique used in the computer program is 
detailed in reference 8. The value of ACL T is found from a first-order truncation of 
the series expansion for Cp in a perturbated flow field which is doubled to account for 
pressure wave reflection from wing surface (refs. 6 and 8). When only parts of the field 
where no shocks are present are considered, ACL T is simply 

ac P,i = j£ r(t) <10) 

and is calculated at the trailing edge of each grid element block within the zone of 
nacelle-wing interference. This interference zone is the area bounded by the intersec- 
tion of the shock surface from the lip of the nacelle (no duct flow spillage being assumed) 
or the nose of the pod (an attached shock being assumed) and the wing camber surface. 
Appreciable interference force error can occur if portions of the interference zone extend 
forward of the leading edge or outward of the wing tip and the disturbance field can spill 
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over onto the upper surface. In such cases, the field point pressure is influenced by both 
the camber surface and the nacelle within the forward-facing characteristic surface, and 
cannot be calculated by the method outlined in this report. 

A mean AC n T over each grid element area is found from 

- c*ep,i (L< - Mq. - 

Outside the interference zone, AC p> i(L*,N*) av = 0. Overall values of interference lift 
and drag are obtained from 


L max N max 

£ £ ACp }I (L*,N*) a v A(L*,N*) 

_ L=1 ^ m in , 


^max ^max 

^ AC p> i(L*,N*) av A(L*,N*)||(L*,N*) 

L.1 N mln 

;3S 

Lmax N m ax 
£ £ B(L*,N*) 

g _ L=1 N m j n 

13 

The weight factors A(L*,N*) and B(L*,N*) are determined from the conditions 


A(L*,N*) = 0 

(l* - x si ^ 0) 

A(L*,N*) = L* - xsi 

(0 < L* - x S i < l) 

A(L*,N*) = 1 

(L* - x s i = l) 

B(L*,N*) = 0 

(L*-x Ze iO; L* - x te = l) 

B(L*,N*) = L*-x Ze 

(0<L*-x^ e <l) 
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B(L*,N*) = 1 - (L* - x te ) 
B(L*,N*) = 1 


(0 < L* - x te < l) 
(l* ~ x iq = i; L*-x te ^o) 


Reflex Camber Surface Determination 


A method for computing a camber surface for a specified loading is outlined in 
reference 1 and is based on the equations 

9z , . i 3 Ar, t \ 1 ff (x - £) AC p (£,? 7 ) d£ dT 7 

|J(x,y) = - J AC p (x,y) + jr P ~ = ( n > 

(y - t?) 2 /(x - 0 2 - /3 2 (y - 1?) 2 


: (x,y) = j* |j(x,y) dx 


( 12 ) 


where ACp(x,y) and ACp(£,??) are the incremental lifting pressures on the wing camber 
surface. These equations are suitably modified to fit the grid-element coordinate system 
shown in figure 3, and adapted for a numerical solution on the digital computer. 

By changing the sign of AC p (x,y) and ACp(|,T?) and using them to represent the 
nacelle interference pressure field, the equations and method can be used to generate a 
reflex surface which will cancel all or a part of the initial interference lift. New force 
characteristics are then determined by computing the interference lift and drag on the 
reflexed camber surface. The required force condition is 

|c L ,c + (1 - K)C L 7| = (C L c + C L>I ) (13) 

— 'unreflexed wing reflexed wing 


that is, the sum of the wing camber lift and the nacelle -wing interference lift minus the 
amount of interference lift to be canceled on the initial wing surface is equal to the 
camber lift plus the interference lift on the reflexed wing surface. A second constraint 
is that the moment change on the unreflexed wing must be equal to that induced on the 
reflexed wing. An iteration method is used to find a solution which satisfies the lift and 
moment constraints. 

The necessity for using both of these requirements can be illustrated by considering 
the 100-percent reflex case. With all the initial interference lift canceled, the C m o 
and the total lift of the nacelle-wing combination at zero angle of attack must be the same 
as those of the unreflexed wing alone at the same attitude so that the initial optimized 
lift and moment loading can be assumed to have been restored. 
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Sonic-Boom Effective-Area Contribution 

Part of the interference lift computation consists of determining what percentage 
has been developed at stations along the X-axis. It is this accumulated lift that deter- 
mines the interference lift contribution to the sonic-boom pressure field directly below 
the airplane. The sonic-boom program requires area distributions which are composed 
of volume and lift components superimposed at their proper locations along the aircraft 
reference line. Interference lift contributes an effective area defined by 

A e (x) = | C L (x) S (14) 

where Cl(x) is the portion of the total interference lift that has accumulated at the 
length x. 


EXAMPLES 

Interference Pressure Determination 

A simplified ducted body (fig. 5) was used to compare numerical and experimental 
values of free stream Cp at a free-stream Mach number of 2.96. The truncated- cone 

forebody had a lip angle ^tan -1 ^ at x = oj of about 3.61° and was 1.75 inches (4.445 cm) 

long. A 2.25-inch (5.715-cm) cylindrical aft section completed the model which although 
a bit stubby, was still a reasonably good nacelle representation. 

The results of wind-tunnel tests and numerical-method calculations are shown for 
comparison in figure 5 for two distances measured in terms of nacelle body length. A 
similar comparison for a slender cone-cylinder body at M = 2.0 is shown in figure 5 of 
reference 6. In both cases, the numerical solution gives a good estimate of the pressure 
signature. The rounding of the nose shock in the ducted nacelle case is due to tunnel 
turbulence, model vibration, probe vibration, and probe boundary layer. 

Figures 6(a) and 6(b) show convergence results for a similar truncated-cone, 
cylindrical afterbody, ducted nacelle located under a rectangular flat plate in M = 2 
flow. For this case which is typical of many design situations, a satisfactory prediction 
is achieved with the use of 20 or more elements along the semispan and 16 or more body 
stations. There is some uncertainty in fixing the minimum number of body stations, 
however, due to the nonuniformity of the wing-intercepted flow field. Both positive and 
negative C p regions occur, and the change in C p across the shoulder is very pro- 
nounced. A sufficient number of F(t) points (corresponding to the number of body 
stations) must be available to define the shoulder expansion and recompression pressures 
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at intervals compatible with the grid element size, or the estimates of lift, drag, and 
moment will have significant errors . 

Interference and Reflex Determination 

A wing-nacelle arrangement (fig. 7) was studied to determine the effect of nacelle 
interference and camber-surface reflexing on performance. The wing was a delta plan- 
form of /3 cot A = 0.5 whose surface was twisted and cambered to provide minimum 
drag due to lift at M = 2.0 for a design Cl of 0.1. Two bodies of revolution composed 
of a truncated- cone forebody and a cylindrical aft section represented the ducted nacelles. 
The forebody had a length of 0.1875?, a minimum diameter of 0.0344?, and a maximum 
diameter of 0.0625?, whereas the cylindrical aft section had a length of 0.094?. 

Figure 8 shows the effect of nacelle interference on the optimum camber surface 
drag polar when no reflexing is present. The addition of nacelles to the optimized camber 
shape resulted in a negligible improvement to the drag characteristics at design Cl> 
but did increase the nosedown pitching moment, as was expected. There may be some 
question concerning the addition of an arbitrary loading (due to nacelle interference) to a 
presumed optimum wing loading to obtain a further drag reduction. However, it must be 
remembered that the nacelle is an external element not included in the original wing 
design problem. In effect, a new optimizing problem is created and an important question 
is raised: In this new situation, will the restoration of the optimum wing loading, or some 
other loading, result in further benefits ? 

Varying the amount of reflexing gave a small but noticeable improvement in drag at 
Cl = 0-1 for relatively small departures from the original optimized camber surface 
(fig. 9). A quantitative analysis to establish a level of reflexing for minimum Cj) at 
design Cl is outside the scope of this report. Figure 9 is included to demonstrate that 
the computer program could be useful in such a study. 

In figure 10, three schematic sketches show, by using the right-hand side only, the 
nacelle-wing configuration, the induced pressure field, and the resultant reflexed camber 
surface in a somewhat exaggerated scale. Although the interference pressures start 
abruptly and change markedly across the chord, the reflex surface varies smoothly and 
gradually, a slope discontinuity occurring only at the border of the interference region. 
The small differences between the original and the reflexed camber can be seen on the 
sketch to the right of the reflexed wing surface drawing. However, the gradual changes 
in slope and height of this reflex surface do not remove the need for the designer to 
exercise discretion and judgment, based on practical experience, in positioning the 
nacelles and selecting the amount of interference lift cancellation. 


11 



The final example (fig. 11) shows how the method is used to provide inputs to a 
sonic -boom analysis of an aircraft configuration. A simple, ducted, parabolic body was 
placed below a flat-plate representation of a wing. The curves described the pressure 
disturbances due to nacelle interference felt directly below the aircraft in terms of an 
effective area ratio. With the nacelle below the wing, the effect is almost double that of 
the nacelle volume taken by itself. Nonlinear behavior of the Whitham flow equations 
causes the ratio of the effective area due to lift to the maximum body area to exceed 1.0. 
The relatively small amount of residual effective area that remains beyond the wing 
trailing edge or the body end point is the result of small errors in the numerical calcula- 
tions and of cutting off the solution at twice the body length instead of at infinity. These 
errors would be inconsequential in the total sonic-boom analysis. 

CONCLUDING REMARKS 

A numerical method for calculating (1) the inviscid flow, nacelle-wing interference 
loads on a warped wing surface, (2) a load- compensating reflex camber surface for opti- 
mization purposes, and (3) the interference-lift effective area contribution to sonic -boom 
pressure disturbances has been presented and discussed. Within the limitations of a 
modified supersonic linearized theory, the numerical method can provide reasonably good 
estimations of the loads imposed on a wing surface by nacelles mounted nearby. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., November 26, 1968, 

126-13-02-08-23. 
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Figure 1.- Pictorial representation of nacelle-wing interference problem. Shaded area indicates interference region. 
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Figure 9.- Effect of reflexing on aerodynamic characteristics of twisted and cambered wing-nacelle combination. 
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Figure 11.- Example of development of nacelle effective area distributions for use in sonic-boom calculations. 
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